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We perform simulations of an effective theory of SU(2) Wilson lines in three dimensions. Our 
action includes a kinetic term, the one-loop perturbative potential for the Wilson line, a non- 
perturbative "fuzzy-bag" contribution and spatial gauge fields. We determine the phase diagram 
of the theory and confirm that, at moderately weak coupling, the non-perturbative term leads to 
eigenvalue repulsion in a finite region above the deconfining phase transition. 

PACS numbers: 12.38.-t, 12.38.Gc, 12.38.Mh 

I. INTRODUCTION 

Recent results from the Relativistic Heavy Ion Collider (RHIC) demonstrate qualitatively new behavior for heavy- 
ion collisions at high energies. RHIC appears to probe a regime where quarks and gluons are no longer bound into 
hadrons but form a new state of matter, a Quark-Gluon Plasma (QGP). The experimental results suggest that the 
transition does not occur directly from a confined phase to an ideal QGP but that it goes through a phase where the 
deconfined phase is non-perturbative In fact, "partial" deconfinement is expected from the Gross- Witten matrix 
model for SU(iV) Wilson lines at large N Q and from quasi-particle approaches Here, we investigate an effective 
theory of straight, thermal Wilson lines [1] which respects the global Z(N) symmetry of Yang-Mills and might be able 
to describe the region below about T ~ 3Td, where resummed perturbation theory fails. We focus, in particular, on 
the distribution of the Eigenvalues of the SU(2) Wilson line at, and slightly above, the deconfining phase transition. 

II. EFFECTIVE THEORY 

We consider an effective theory in three dimensions which uses the thermal Wilson line L(a;) as the fundamental 
degree of freedom in the electric sector. This is a matrix valued field which exists in three dimensional space. The 
effective Lagrangian is composed of a kinetic term, the one-loop potential, and a non-perturbative term proportional 
to a "fuzzy-bag" constant Bf. 



£ cff = ~tr G% + |Lt A:L| 2 - 4,T 4 V ^|trl/f + B/T 2 |trL| 

n>l 



Gij is the non-Abelian magnetic field strength and Dj the covariant derivative. This theory does not involve elements 
of the algebra, i.e. A , and respects the Z(N) center symmetry of the SU(N) gauge group. It is therefore able to generate 
fluctuations between different Z(N) vacua, just above the phase transition. The Lagrangian is non-renormalizable in 
three dimensions and thus applies only at distances 3> 1/T. 



III. LATTICE SIMULATION 



We performed simulations of the effective theory for two colors, employing three-dimensional cubic lattices with 
periodic boundary conditions. We use Metropolis and overrelaxation techniques [6, 7] to update the lattice configura- 
tions and to generate a thermal ensemble. We currently neglect all but the n — 1 term of the perturbative potential 
which can then be combined with the non-perturbative term. The lattice action then reads 

S = 0Y,( 1 ~ t^RcTi'Uq) - i/?£> (LiU«L}uJ + h.c.) - m 2 ^ |trLi| 2 . 

□ <ij> i 
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FIG. 1: Eigenvalue distributions measured on a 24 3 lattice for different values of m 2 and (3. Magnetic fields are included. 



It contains the standard kinetic Wilson action for the magnetic fields, where the sum runs over all " plaquettes" , 
gauge-covariant kinetic nearest neighbour interaction for the Wilson lines and a mass term. The index i labels lattice 
sites and (ij) labels links. 

Initial results were obtained by neglecting the magnetic sector [5[- The theory is then essentially reduced to a 
sigma model with global SUl(2) x SU#(2) symmetry. When m 2 = 0, a second order phase transition occurs at 

[3 C = 0.942(5) . Here, the expectation value of the "length" of L, u = \J tr L^L/2, takes on a non-zero value 8]. The 
effective screening mass, which equals the inverse spatial correlation length, drops to zero here (up to finite lattice 
size effects) and remains zero above the phase transition point, due to massless Goldstone modes. 

Choosing a fixed value for (3 and varying m 2 also gives a second order phase transition. Here, the proper order 
parameter is the expectation value of the Polyakov loop, which is the trace of the Wilson line. Larger (3 shifts the 
phase transition point to smaller m . To obtain non-zero VEVs for the Polyakov loop for negative m 2 , an infinitesimal 
background field must be added. These results suggest that the phase boundary of the theory runs diagonally in the 
(3-m 2 plane 

We measured the average and the difference of eigenvalues of the Wilson lines 



pi(t,x) = ^ \Xi(t,x) - X 2 (t,x)\ 



P2(t,x) = i \Xi(t,x) + A 2 (t,x)| 



and determined their distribution in the thermal ensemble over a broad range of couplings. We find, that for m 2 = 
and small /3, the distribution is dominated by the group integration measure which leads to a logarithmic divergence 
for pi — * 0. For large m 2 one approaches the perturbative vacuum, where the distribution of the average p 2 is peaked 
around 1. In the confined phase at large (3 the Wilson lines fluctuate about the non-trivial vacuum L c = ir^, or SU(2) 
rotations thereof (eigenvalues "repel" ). Fluctuations diminish with increasing /3. 

Adding gauge fields now leads to a few qualitative changes in the observables. Foremost, it is no longer necessary 
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FIG. 2: The spatial string tension as a function of (3 (left) and as a function of m 2 for fixed /3 = 2.0 (right). The string tension 
appears to be independent of m 2 . 



to use the 0(4)-length as an order parameter, since the presence of gauge fields explicitly breaks the SUl(2) x SU,r(2) 
symmetry of the kinetic term. The Polyakov loop is now a good order parameter, even for m 2 = 0. All massless 
Goldstone modes disappear as well, so that the inverse correlation length is only zero (again up to finite size effects) 
exactly at the phase transition point. 

Qualitative differences in the distribution of eigenvalues arise from the fact that even at m 2 = raising (3 can now 
lead to eigenvalue attraction. The qualitative dependence on m 2 for fixed /3 is similar to the spin-model case, although 
the transition from the perturbative vacuum to the confined vacuum is smoother across the point m 2 = (as one 
would expect, since this is reminiscent of the change of the behavior of the Polyakov loop). Overall, changes in the 
electric sector due to magnetic fields are qualitatively similar to the effect of the external field which was formerly 
used to find the phase boundary in the spin model at m < 0. 

To assess the influence of the dynamics of the Wilson lines on the magnetic fields we have also measured the 
expectation value of the Wilson action itself and of the spatial string tension. This was done by obtaining expectation 
values of rectangular Wilson loops W(I, J), where I and J are the dimensions of the loop in lattice units. Such loops 
obey the area law W(i,j) oc exp (—a A), where A = a 2 IJ is the area of the loop and the string tension a plays the 
role of a decay constant. One normally uses Creutz ratios [9( but our comparison of both methods suggests that in 
3D, using the Wilson loops directly works equally well. 

We find that both the string tension and the Wilson action drop approximately exponentially with /3 but are 
independent of m 2 . Both are also independent of the lattice size. The magnetic sector is largely unaffected by 
deconfinement in the electric sector. 
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